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ABSTRACT 

Studying the environments of 0.4<z<1.2 UV-selected galaxies, as examples of 
extreme star-forming galaxies (with star formation rates in the range of 3— SOM© 
yr^^), we explore the relationship between high rates of star-formation, host halo 
mass and pair fractions. We study the large-scale and small-scale environments of 
local Ultraviolet Luminous Galaxies (UVLGs) by measuring angular correlation 
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functions. We cross-correlate these systems with other galaxy samples: a volume- 
limited sample (ALL), a Blue Luminous Galaxy sample (BLG) and a Luminous 
Red Galaxy sample (LRG). We determine the UVLG comoving correlation length 
to be ro = 4.8^2 4^ h~^Mpc at (z) =1.0, which is unable to constrain the halo mass 
for this sample. However, we find that UVLGs form close (separation < 30 kpc) 
pairs with the ALL sample, but do not frequently form pairs with LRGs. A rare 
subset of UVLGs, those with the highest FUV surface brightnesses, are believed 
to be local analogs of high redshift Lyman Break Galaxies (LBGs) and are called 
Lyman Break Analogs (LBAs). LBGs and LBAs share similar characteristics 
(i.e., color, size, surface brightness, specific star formation rates, metallicities, and 
dust content). Recent HST images of z~0.2 LBAs show disturbed morphologies, 
signs of mergers and interactions. UVLGs may be infiuenced by interactions 
with other galaxies and we discuss this result in terms of other high star-forming, 
merging systems. 



Introduction 



The observed properties of galaxies can be attributed to a wide variety of complex 
factors. In order to understand the environmental causes that might trigger intense star- 
formation in the local universe, such as that observed at high redshifts for Lyman Break 
Galaxies, we explore two regimes: large-scale, to study the connection between halo mass 
and high star formation rates (SFRs); and small scale, calculating merger rates. At large 
scales, a galaxy's SFR may be connected to how the galaxy accretes gas from its large-scale 
environment, which depends on the halo mass. Smaller scale environmental effects, from 
nearby galaxies, may induce star formation from interactions and eventual mergers. 

The study of galaxy evolution has been connected to that of larger structure formation- 
their growth related to overdensities arising from primordial density fiuctuations. Studies 
of galaxy bias, the distribution of galaxies relative to underlying dark matter distribution, 
attempt to describe the fate and origin of present-day galaxies. The link between halo mass 



and star formation has been discussed and explored in several papers (IWang et al.l 12006 



Mandelbaum et al.l l2006l : IConroy et al.l l2008l : iKhochfar fc Ostrikerl |2008| ) . A simple descrip- 
tion involves gas flowing to centers of dark matter potentials, reaching sufficient density to 
radiate and cool, losing pressure support and settling into a disk, and then forming stars. 
The system is eventually accreted onto a larger structure. Such accretion by clusters can 
shock-heat the gas into high temperatures, preventing new stars from forming. Therefore, 
we expect that active star formation can occur for halo masses whose virial temperatures 
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are less than lO^K, where shock-heating inhibits new star formation. 



Redshift evolution of the link between halo mass and star formation, determined by 
l arge scale en vironmental effects, might explain observations such as the Butcher-Oemler 
(jOemlerlll974l ) effect. In this effect, clusters at high redshifts host more blue, star-forming 
galaxies than local clusters, which are predominantly inhabited by "red and dead" galaxies. 
Similarly, us i ng th e two-point angular correlation function to measure clustering strength, 
Heinis et al.l (120071 ) find that UV-selected, star-forming galaxies have migrated from cluster 
regions, with high correlation lengths (tq), at high redshift to more isolated regions at low 
redshifts. From the view of hierarchical structure formation, higher redshift clusters have 
lower halo masses compared to low-z clusters, thus high-z cluste rs have viria l temp eratures 
low enough to allow star formation. Adding to this framework, iKeres et al.l (120051 ) discuss 
the importance of hot vs. cold accretion - in the simphfied picture above, gas does not get 
heated to the virial temperature necessarily. Massive galaxies draw gas via quasi-spherical, 
"hot mode" accretion, in which the gas does shock-heat to the virial temperature. However, 
less massive galaxies acquire gas from flow along filaments in "cold mode" accretion where 
the gas radi a tes away its gravitat i onal e nergy and never reaches the virial temperature (see 
Keres et al.l (120051 ): ICroton et al.l (120061 ) and references therein). These simulations show 
that the large-scale galaxy environment influences gas accretion processes. Therefore, a 
galaxy's halo mass contributes significantly to how it acquires gas, the fuel that affects its 
star formation efficiency. 

Similar sho ck physics may affect the relation between small scale merger effects and 
star formation. iKauffmann. White, fc Guiderdonil (119931 ) predict that the merger of a halo 
with a larger structure, such as a cluster, suppresses star formation as gas is shock-heated 
to the virial temperature of the more massive halo and then cools onto the central object of 
the most massive halo. Thus, these results suggest that the mass of the larger companion 
is a co nsideration for whether mergers raight e nhance star formation. According to the 
work of IKauffmann. White. &: Guiderdonil (Il993l ). mergers with clusters should not result in 
enhanced star formation. 



Nevertheless, other studies show that at these small scales galaxy-galaxy interactions 
and close neighbors can stimulate star formation activity. In hierarchical structure formation, 
the merger process is an important ingredient for galaxy evolution. Over thirty years ago, 
Toomre fc Toomrd (119721 ) suggested that tidal forces during the merger process can funnel 
gas into the central regions of galaxies, triggering strong star formation activity. Further 
theoretica l studies have invest i gated t he connection bet ween mergers and high star fo rmation 



rates (see 



Mihos fc HernquistI (Il9941 ). ICox et al.l (l2006l ). and ldi Matteo et al.l (120071 )) 



Li et al.l (120071 ) have shown that nearly half of the highest star-forming galaxies have 
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close companions (galaxies within a projected radius of 100 kpc), by correlating 10^ galaxies 
from the Sloan Digital Sky Survey (SDSS). According to this research, tidal interactions are 
the primary cause for enhanced star formation for these most active star-forming galaxies. 
Several studies using various star formation indicators from Ha, CO, far-infrared luminosi- 
ties for interactin g gal a xies have shown tha t enhanced star formation results from galaxy 
interactions (e.g.. iKeell (Il99ll ). IStruckl (119991 ) . and references therein). 



Using the Galaxy Evolution Explorer {GALEX ), iHeckman et al.l (120051 ) have uncov- 
ered a population of local, intensely star-forming galaxies. These galaxies are selected based 
on their FUV luminosities (Lpuv > 2x10^^Lq). The most compact of these UV-luminous 
galaxies (UVLGs), with Ipuv > 10^ kpc~^, share many similar proper ties with z^3 LBGs : 
specific SFRs, metal l icities, morphologies, ki nematics, and attenuations (IHoopes et al.ll2007l : 
Overzier et al.l l2008l : iBasu-Zych et al.l 120071 ) . We refer to these systems as Lyman Break 
Analogs (LB As), to discriminate this special subset of UVLGs from lower surface brightness 
UVLGs. Given that LBGs are found at higher redshifts (z=2-3) where major and minor 
mergers are generally more prevalent than for z=0.1 LBAs, LBAs and LBGs may be ex- 
pected to differ in their recent merger histories, environments, and location within large 
scale structure. 

In this paper, we study the environments, as quantified by correlation functions, of 
low to intermediate redshift (0.4<z<1.2) UVLGs in the COSMOS field. We extend our 
analysis to LBAs to compare with z~3 LBGs, although this analysis is limited because of 
the rarity of z<l LBAs. We measure two-point auto-correlation functions (ACFs) and cross- 
correlation functions (CCFs) of UVLGs with other galaxy samples to probe the mass of the 
dark halos within which UVLGs reside. We connect these observations with those of LBGs 
at high redshift (z~3), which are known to strongly cluster ( Adelberger et allboosl ). Recent 
Hubble Space Telescope (H ST ) images of z~0.2 LBAs show disturbed morphologies, signs 
of mergers and interactions (jOverzier et al.ll2008l ). We relate small-scale correlations to pair 
fractions in order to estimate merger rates. Studying the locations of current UVLGs, as 
examples of extreme star-forming galaxies, we explore the relationship between high rates of 
star- formation, host halo mass, and merger rate as a function of redshift. 

In the following section, we discuss our sample selection. In the three subsequent sec- 
tions, we outline the theory ( §2.ip and techniques for measuring the angular correlation 
functions( §2.2p . and derivations of physical quantities ( §2.31) used to analyze our data. We 
discuss our results for the large and small scale environments of UVLGs and LBAs in ^ 
and discuss these in the context of other environmental studies in ^ A summary of our 
analysis can be found in the final section: ^ Throughout this paper we adopt the cosmology 
given by (Hq, ^m, ^a) = (70 km s"^ Mpc~^, 0.3, 0.7). Unless explicitly stated, our number 
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densities are in physical, proper units. Correlation lengths are quoted in comoving units of 
h-^Mpc. 



2. Data and Sample Selection 



With its large contiguous field and deep observations, the COSMOS survey provides 
excellent data for environmental studies. The COSMOS field covers 2 deg^, centered on 
aj2ooo = 10'*00"'28.6^ 5j2ooo = +02°12'21.0". This region of sky has rich multiwavelength 
data, including: optical observations with the superb Subaru ground-based telescope, atop 
Mauna Kea, Hawaii, ultraviolet coverage from GALEX , radio observations from VLA, and 
Xray data from XMM. The inner 1.7 deg^ has been imaged by HST , u sing the Advanced 
Camera Survey (ACS) in the F814W filter (see iKoekemoer et al.l (120071 ) for details). Deep 
photometric data from Subaru provides ph otometric redshift s with Az/(1 + z) < 0.02, de- 
pending on I-band magnitude and redshift (lllbert et al.ll2008l ). 



From the I-band selected (I<25 . 5) CO SMOS photometric redshift catalog (for details 
about this catalog see ICapak et al.l (120071 )). we compare UVLG and LBA s amples, which 
draws from the combined GALEX data (I<25.5, and NUV<23, See lZamoiski et al 

(120071 ) for more details), with three separate galaxy samples (selected from the COSMOS 
sample, with I<25.5) for 0.4<z<1.2: (1) a volume limited sample in Mg-"ALL", (2) Blue 
Luminous Galaxies- "BLG" , and (3) Luminous Red Galaxies- "LRG" . We briefiy note here - 
while the "ALL" sample is volume-limited in Mg, it is not volume-limited for all colors. This 
is discussed in more detail further below. Figure 1 shows the various cuts applied to the data 
to create these separate samples: ALL is comprised of galaxies with -24< Mg <-18; BLG 
galaxies have Mg <-20.0 and M^ — Mr < 1.8 color; and LRG set is selected by a luminosity 
cut, Mr <-20.5 and color cut, Mu — Mr > 2.0. Both colors and luminosities are given in rest- 
frame, k-corrected AB mag nitudes. K-corrections were done using the kcorrecL4-l-4 IDL 
routine (IBlanton et al.ll2003l ). The spectral energy distribution (SED) fitting was done using 
the COSMOS optical bands as well as the GALEX FUV and NUV bands, where available- 
since the comparison samples are not required to be UV-detected. We modified this routine 
to use the closest band to the rest-frame filters for determining kcorrections- for example, 
at z=0.5, the observed NUV magnitude at 2300A best approximates the rest-frame FUV 
magnitude at 1500A(i.e, 2300A = 1500A * (1 + 0.5)). 

The UVLGs are selected using their k-corrected, rest-frame FUV properties from the 
COSMOS+ GALE X catalog: FUV <-1 9.54, corresponding to the luminosity criterion of 



Lpuv > 10 Lq (iHeckman et al.l 120051 ). To estimate the FUV surface brightness, we use 



the value of rso output from SExtractor (IBertin fc Arnoutslll996l ) for the ACS I-band data 
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(see IZamojski et al.l (120071 ) for further detail), which corresponds to the half-hght radius. 
We use this ACS I-band half-hght radius as a proxy for the UV size, acknowledging that 
the I-band at these redshifts corresponds to optical wavelengths. Our sample of LBAs were 
determined by the criteria: Ipf/y > 10^ Lp kpc~ ^. While the I-band size may overestimate 



the UV size by a factor of 2 (jOverzier et al.l (120081 ). thereby, underestimating the FUV surface 
brightness by a fac tor of 4), we find th at ~ 20% of UVLGs are LBAs, consistent with what 
is found at z~ 0.2 (IHoopes et al.l 120071 ). 



We summarize the sample sizes for each redshift bin in Table 1, and show the redshift 
and spatial distributions of these samples in Figures 2 and 3, respectively. 

The "L" shape of the color-magnitude diagrams in Figure 1 indicates that the red and 
blue galaxies do not populate the magnitude space in the same way. The majority of red 
galaxies tend to be more luminous than blue galaxies; the faint end for the red galaxies starts 
decreasing around Mg =-19. While the volume-limited cut shown on the upper left panel of 
Figure 1 selects a complete sample in Mg, the color-magnitude diagram at Mg ~-19 shows a 
diminishing sample of red galaxies. The lur ninosity funct i on (L F ) for red galaxies sh ows that 
the faint end does have a decreasing slope (IZucca et al.l (120061 ): iFaber et al.l (120071 ) . but our 
LF comparison of the red galaxies in the ALL sample with the theoretical LF does indicate 
sample incompleteness at higher redshifts. This incompleteness does not affect either the 
BLG or LRG samples which are selected to be sufficiently bright to avoid magnitude limits. 
However, the ALL magnitude cut of Mg =-18 suffers from some incompleteness in the red 
population. This observation is worth noting for interpreting and comparing the ALL sample 
to other samples and studies, but does not affect the overall conclusions of this paper. 

The predicted n(z), shown as a black dotted line in Figure [21 is derived using the 
lumino sity functions observe d by iFaber et al.l (120071 ) for the ALL, BLG and LRG samples 
and by lArnouts et al.l (120051 ) for the UVLG sample (this is discussed further in §2.3p . The 
redshift distributions in Figure 2 roughly follow the smooth predicted distribution marked 
by the black dotted line. The UVLG, as well as the LBA, redshift distributions may be 
too noisy, because of poor number statistics, to interpret their shape. The overdensities 
apparent in the spatial maps (Figure Bj) corr e spond to known physical s tructures obse r ved in 
the COSMOS field (see IScoville et air(l2007h , iFinoguenov et al.l ( l2007h , iMassev et all J2007h 
for more discussion of these structures). 
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2.1. Correlation Functions — the Theory 

In the standard framework for defining the two-point correlation function, galaxies are 
distributed in space such that the probability, 6P, of finding a galaxy in some volume, 6V is, 
6P = n6V, where n is the mean number density of galaxies. Extrapolating to the situation of 
finding two galaxies at some comoving separation, ri2, one within volume, 6Vi and another 
within volume, (5V2, is then 

6F = n^5Vi5V2[l + ^iTu)] (1) 

where ^(ri2) is defined as the two-point correlation function. An approximation for the 
correlation function, ^(r), is a power-law with a single slope at all scales: 

e(r,z) = (r/ro(z))-^ (2) 

where z refers to redshifft, 7 and ro(z) specify the slope and comoving clustering strength, 
respectively. Given a uniform, Poissonian distribution, ,^(r) = 0, while correlated data 
(^(r) > 0) means that there is an increased probability of finding a second galaxy within 
distance r of the first. 



Various studies (i.e.. IZehavi et al.l (120041 ) for SDSS galaxies and lOuchi et al.l (120051 ) for 



high-z LBGs) have found departures from a power-law with a single slope. Rather, they find 
that the correlation function is better approximated by the sum of two contributions: the 
"one-halo term" , dominating at small scales for galaxy pairs residing within the same halo, 
and the "two-halo term" , dominating at large scales for galaxy pairs within separate haloes. 
This so-called HOD modeling is beyond the scope of this paper. However, to account for the 
fact that the galaxy correlation function might not follow a power law with a single slope, 
we assume that ^(r) can be fit by two power laws with different slope for small scales and 
for large scales. 

For data based on photometric redshifts, the three dimensional correlation function can 
not be measured directly, but it can be inferred from the projected two-dimensional angular 
correlation function. By analogy of the steps outlined above, the three-dimensional, real 
space ^(r) is collapsed into a two-dimensional function: w(0), which measures the excess 
probability that two galaxies are found within some angular separation, 6. 

5F = N'^[l + w{e)]5ni5n2 (3) 

where N is the surface density and 6Qi and 6Q2 are the solid angle elements describing our 
field. Like ^(r) in Eqn. [21 w(^^) is also expected to have a similar power law form: 

w(^^) = a^,eO-'. (4) 
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We convert the three- dimensi o nal co rrela tion fu n ction into the angular one, using the for- 
malism introduced by iLimberl (119541 ) and iRubinI (jl954l ). Therefore, we relate the comoving 
correlation length, ro, to a„^e and 7 to 5 by the following: 



7 

g(z) 



5 + 1 

/^ro(z)^g(z)(^i»2)dz 
Xni(z)dz/^n2(z)dz 



(5) 
(6) 

(7) 



where = Jj^^dx(l + x^) '''/^ = ^^^7)^ \ r is the comoving distance, ni{z) and n2{z) 
are the number densities, and F(r) is the curvature term from the Robert son- Walker metric, 

ds^ = cV - a2[drVF(r)2 + r^dO^ + r^sin^^-^-^^i 



Using the iGroth fc Peebles! (119771 ) form for the evolution of ro(z) 

ro(z)=ro(0)(l + z)-(3— 7)/7 



(9) 



where ro(0) is the comoving correlation length at z=0, and e gives the evolution of clustering 
(i.e., e = specifies a stable state, e = 7 — 3 describes no evolution, and e = 7 — 1 corresponds 
to the linear growth of fluctuations). Throughout our analysis, we assume no evolution, or 
e = 7 — 3, which allows us to simplify Eqn. O Over the relativ ely narrow redshift ra nges 
for which we are integrating over, we can justify this assumption- iMeneux et al.l (120061 ) find 
that ro evolves by < 20% between z=0.2 and z=1.2 for late type galaxies and < 50% for 
early type galaxies. 

Auto-correlations compare distributions of a sample with itself, while cross-correlations 
compare distributions of one set of objects with a separate sample. Naturally, in the auto- 
correlation cases, ni = n2 in Eqn. |6l We infer the real-space correlation length, ro and 
power-law slope, 7, by fitting the angular correlation function. 

Along with slope and amplitude parameters, we additionally calculate the Integral Con- 
straint (IC) correction. In principle, our observed field may be superimposed on a void or 
overdensity, such that the observed number density in our field will not match the average 
number density in the universe. In theory, the integral in Eqn. [3] should average to over all 
scales, i.e., since there are N galaxies in the universe, excess at small scales must be balanced 
by deficits at lar ger scales. We c orrect the observed angular correlation function by the IC 
factor defined in lRoche fc Ealed (119991 ). 



- 9 - 



2.2. Measuring the angular correlation function 

The angular correlation function, w(^), measures the excess of clustering compared to 
a uniform, random distribution. To simulate the random distribution appropriately, we first 
create masks from the COSMOS+ GALEX fields. These masks mark the locations where 
bright, saturated stars and field edges prevent galaxies in those locations from appearing 
in the catalog. The black regions in Fig. [3] display masked locations. Populating the 
allowed, un- masked regions with random galaxies, we create the Random set. The histogram 
of separations in the Random set and Data sets are referred to as RR(0) and DiD2(^), 
respectively, where the subscripts 1 and 2 denote two separate data samples; D1D2 =DD in 
the auto-correlation case. Similarly, the separations between galaxies in the Random and 
Data sets are represented by DiR and D2R, or simp ly DR for autocorrelations. It has been 



shown that the Landy-Szalay (iLandy fc Szalaylll993l ) estimator is one of the most robust ways 



of measuring w(^), reducing the edge effects of finite fields. Here we give the expression for 
the angular correlation function that we use to measure auto-correlations (Eqn. [TUj) and 
cross-correlations (Eqn. [TT]) : 

^ ^ n(n - 1) RR n RR ^ ' 

r(r-l)^DiD2 ui D2R n2 DiR ^ nin2 RR\ ^^^^ 



nin2 \ RR r RR r RR r(r - 1) RR 

where r and n refer to the number of galaxies in the Random and Data samples (the sub- 
scripts 1 and 2 for ui and n2 distinguish the separate samples, in the cross-correlation case). 
O ur angular correlation function error s have been derived from boot-strapping, as described 



m 



Barrow. Bhavsar. fc Sonodal (119841 ). Cosmic variance adds another important source for 



error. However, unlike the bootstrap errors which ar e scale dependent, the errors from cos mic 



variance shift the entire correlation function. Using ISomerville. Primack. &: Faben (120011 ) to 
calculate cosmic variance errors for our samples, we show these errors in our correlation 
function figures as shaded regions (for example, see Figs. [Mil)- 

As described in ^ our samples span the redshift range 0.4<z<1.2. While we may 
reduce statistical noise by choosing the largest sample possible - selecting galaxies in this 
entire redshift span, correlating data across such a wide redshift range has disadvantages. By 
correlating data in two- dimensions, which aren't physically correlated in three-dimensional, 
de-projected space, we risk weakening the actual correlations. Therefore, we select two bin 
sizes: dz=0.2 (resulting in 4 overlapping bins: 0.4—0.6, 0.6—0.8, . . . , 1.0—1.2) and dz=0.4 
(2 bins: 0.4—0.8, and 0.8—1.2). Comparing the angular correlation functions for different 
redshift bins requires converting the angular correlation function into physical de-projected 
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units, by the following reasoning: 



dN(n) /dN^ 



, i) 

= dN^"^° (12) 

\dn/o 



+00 

21 1/2 



00 



dN(rp) _ /dN\ 
dA \dA/0 ^-[^^-j 



\dV/o 

where rp is the transverse distance and n is the line-of-sight separation; (^)q, (^)o; 
(sv)o I'sspectively, dN/dA, dN/dfi and comoving number density for a uniform, non- 
clustered distribution. Finally, multiplying the w(^) equation (Eqn. [T2l) by dQ/dA and 
averaging these expressions over the redshift bins, we derive 



w(^) _ (f> 



Wn r- 



/ dN _ dn\ 



p^"p^ \dn dA/ 

dn 



dA D2(l + z)2 



(15) 
(16) 



dN 1 

Wp(rp) = t^nT w(e) (17) 

\dV/ 

where Da is the angular diameter distance. As such, we can compare how the correlation 
function of the narrower bin sizes compare to the wider redshift ranges. Fig. H] shows how 
the wide bins differ from the narrow ones for samples with 0.8<z<1.2; we plot the ratio of 
the correlation functions for the wide to narrow bins, displaying both w(^^), top, and Wp(rp), 
bottom. Since the discrepancy was not severe, we decided that the wider redshift bins would 
provide better statistics. Given the low abundance of UVLGs, their auto-correlation only 
proved possible for the full 0.4<z<1.2 sample; and for the LBAs, their auto-correlation 
proved too noisy to be useful, and we were only able to attain reasonable results by cross- 
correlating with the largest sample: the ALL sample with dz=0.4 bins. 

For each sample, we calculate separate fits to the angular correlation function in two 
regimes: large scales (> 200kpc) and small scales (< 200kpc). We create a grid of a^, spaced 
logarithmically, and delta, linearly spaced, values, with the best fit values corresponding to 
the minimum x^- The errors on this fit are derived as the values corresponding to min(x2) + l, 
and we make certain that the grid values extend significantly past these error values to 
avoid restricting the errors within the minimum or maximum value. The large scale is fit, 
additionally to the slope and amplitude, with the IC correction, discussed in the previous 
section. We show examples of our fitting to the separate samples in Fig. O 
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2.3. Calculating Physical Quantities 



The number density is required for calculating both Limber's Equation and the pair frac- 
tion. For our analysis, we use the observed number density. However, as a check we calculate 
the predicted physical number densities for the ALL, BLG, and LRG samples. We u s e a, M p 
(given in Vega magnitudes, and converted to AB magnitudes using IWillmer et al.l (120061 ) ). 
and 0* given for the COMBO-17 data in iFaber et al.l (120071 ) to determine the luminosity 
function, interpolating the LP between the specified redshift-dependent values to derive a 
smooth n(z) across our entire redshift range. We use a Gaussian kernel with a = 0.02, 
simulating the photo-z errors, to convolve the derived n(z)- the results are shown as black 
dash-dotted lines in Fig. [21 compared to the observed numbers. With the exception of some 
peaks and valleys, the observed BLG and LRG samples follow the predicted values well. 
The agreement between the predicted and observed ALL distribution is poorer- at higher 
redshifts the observed number density appears to be underestimated. It is possible that the 
peaks in the observed n(z) in the LRG sample correspond to real structure. For the UVLG 
and LB A samples, we pre dict the number density using the calibrated luminosity function 
from lArnouts et al.l (120051 ). In Table [21 we list the values used to fit the luminosity function 
to derive n(z) for all of our samples. For use in Limber's Equation, we smooth the observed 
redshift distribution, performing a boxcar average with a kernel size: Az=0.08. The dis- 
crepancy between results derived using the observed versus predicted number densities is 
negligible. 

Quantifying the smaller scales, we relate the fits from the inner regi ons (shown in red 
in Fig. [5l) to pair fractions and merger rate, similar to analysis done by [Bell et al.l (120061 ) 
for massive, luminous galaxies. The pair fraction is derived from integrating Eqn. [Tjfrom 
to some maximum radius, within which galaxies are assumed to eventually merge- the pair 
fraction quantifies the number of galaxies that are found within some critical, merging radius 
of another galaxy. However, since we are interested in the physical separation between the 
galaxies rather than the comoving distance, we replace the co moving correlation length, ro, 
with the physical correlation length, ro/(l + z). Therefore, we calculate pair fractions by 
simp lying the integral of Eqn. [T]with the following assumption: ^(r < rf) >> 1, into this 
expression 

4vrn2 / ro 



P(r < rf) 



i3-7 



3-7 \l + z^ 

Here, i s the number density of the second sample in the cross-correlation. Following 
Bell et al.l (l2006l ). we also use rf=30 kpc to define close, physical pairs. Our errors for 
the pair fraction include cosmic var iance uncertainty, estimated using the relations from 
Somerville. Primack. fc Faberl (l200ll ). While referre d to as a probability , the pair fraction 
is the mean number of galaxies found within rf (See iPatton et al.l (120001 )) and may exceed 
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unity. The pair fraction has quantitative power in comparing small scale clustering between 
different samples. 

The merger rate is derived by dividing the pair fraction by the timescale for this merging 
event to occur (assuming that the close pair merge into one merger remnant, such that the 
pair fraction is equal to twice the merger fraction). Though the calculation is uncomplicated, 
the steps for deriving the na erger timescale are quite complex and fraught with uncertainties. 
Kitzbichler fc White! (120081 ) have compared pair fractions derived from their mock catalogs 
with merger rates from the Millenium Simulations to determine merger timescales. From 
these results, we infer that the merger timescale is ~ 3000 Myr h~^, using the most general 
cross-correlation case of UVLGs with ALL and estimating < M^, >= 10^° Mq (IHoopes et al. 



20071 ) in the range z=0.4— 1.2. Using this value, we d eterra ine merger rates for our data, 
shown in Table HI According to iKitzbichler fc Whitd (120081), merger t i mesca les have been 
underestimated by ~ 10 in previous studies- for example, iBell et al.l (120061 ) estimate the 
merger timescale to be 400 Myr. To avoid the large uncertainties in the merger timescale, 
we compare our pair fraction results to other studies, rather than comparing merger rates. 



3. Results 



Given the fits shown in Fig. [HI we derive the correlation strength (ro) for UVLGs. We 
determi ne that UVLGs have tq = 4.8^*^2^4^ h~^Mpc (see Table |3]), corresponding to Mhaio > 



10^^ (see lMo fc Whitd (120021 )). Figure [8] shows how our results compare with other published 
data, at variou s redshifts. The shaded r egion in this Figure corresponds to the shaded arrow 
in Figure 11 of lAdelberger et al.l (120051 )- this marks how tq evolves with redshift for a halo 
of mass IO^^'^Mq (bottom of shaded region) to lO^^'^M© (top of shaded region). We discuss 
the evolution of UV-selected galaxies, in relation to other samples, in the following section. 

Studying the autocorrelations of the separate samples with each other, we find that 
typically (at all z) LRGs cluster most strongly with themselves at small scales, compared 
to auto-correlations of BLGs or ALL samples (top panel of Fig. [7]). Comparing auto- 
correlations of the samples with cross-correlations with UVLGs in the bottom panels of 
Fig. m we find that at small scales (r^ < 200kpc), UVLGs clu ster less strongly with LRGs 
compared to LRGs with themselves {i.e. iMasjedi et al.l (120061 )): however, it is difficult to 
differentiate the UVLG cross correlations with the ALL and BLG samples compared to their 
auto-correlations. 



We compare the pair fractions of the autocorrelation of the other three samples with 
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their cross-correlation with UVLGs in Fig. [Qr I . At the smallest scales the ALL sample 
pairs more with UVLGs than with themselves, and given the large errors, the BLG sample 
galaxies appear to pair with UVLGs with equal likelihood to pairing with themselves. Yet, 
UVLGs are not likely companions of LRGs. Finally, the LBA-AU (purple) pair fraction in 
Table H] has excessively large errors, preventing a meaningful comparison of the LBA-ALL 
sample with these other cases. 

As we have discussed, small scale clustering translates to merger rates- suggesting that 
UVLGs appear to merge often with ALL sample galaxies (and possibly, with other blue, 
star- forming galaxies), and do not merge often with LRGs. However, as there are some 
small scale pairs of UVLG-LRGs, these systems do interact on some occasions. Given that 
blue, star-forming galaxies typically do not reside in clusters, and UVLGs merge with these 
galaxies as often as they do with themselves, we expect that UVLGs will inhabit lower mass 
halos rather than massive clusters- further constraining our results. 

Comparing with other results, we show our pair fraction re s ults fo r the ALL, BLG, 
and LRG samples with other studies in Fig. [TUJ IConselice et al.l (120031 ) observe a sample 
of Mb < —18 galaxies (orange squares) to measure merger fractions, as sumed to b e hal f 
of the pair fractions since two pair galaxies interact in a single merger; iLin et al.l (120041 ) 
study pair fractions for galaxies with —22 <M r < —20 (green triangles), and the magenta 

(120061 ). Since all of these samples 



star marks the result for Mb < —20, found by iBell et al 



are magnitude limited and color-selection independent, t he ALL sample i s most cornparabl e 
(Mg < —1 8 ). It should be noted that many studies (e.g. iLin et al.l (120041 ). IConselicd (120061 ). 
Bell et al.l ( l2006l )) find that pair fractions are extremely sensitive to the limiting magnitude. 
In order to comp are our pair fraction results with other work, we adopt the observation from 
Bell et al.l ( l2006l ) that the correlation parameters are less sensitive to limiting depth than 
the number density. Therefore, we modify the derived pair fraction for our ALL sample, 
using Eqn. [TS] with the small-scale correlation function parameters in Table El but changing 
the number density to represent Mb < —20 selection. The modified pair fraction is shown 
as open circles in Fig. [10] and agree well with published results. 



^''While the pair fraction approaches and exceeds unity for the ALL auto- and cross-correlation results, 
this does not suggest that ~100% of the galaxies are found in pairs. Rather the pair fraction represents the 
mean number of secondary galaxies likely found near the primary ALL galaxy. See discussion in §2.31 
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4. Discussion 

The picture that has emerged from this analysis depicts UVLGs as galaxies that inhabit 
halos with Mhaio > IO^^Mq. UVLGs are more likely to interact and eventually merge with 
other galaxies (ALL) than the ALL sample does with itself, and they are as likely to interact 
and have companions as other blue, star-forming galaxies. These results are consistent with 
other studies showing that high SFRs may be closely related to increased merger rates. 
UVLGs are unlikely to interact with LRGs, which are expected to reside at the centers 
of massive clusters, providing another constraint regarding UVLG halo masses. In this 
section, we first place UVLGs amongst other categories of star-forming galaxies; then, we 
discuss our findings in light of other work related to UV-selected galaxies, star-formation 
and environments. 



4.1. UVLGs vs. other star- forming galaxies 



In Figure [H we display the large scale en vironmental measure fo r several other classes 
of star- forming galaxies. The BM, BX, LBGs (lAdelberger et al.ll2005l ) are selected by color 
and limiting magnitude to be both luminous and blue, exhibiting the Lyman Break in their 
spectra. As discussed in the introduction, the LBAs have been found to have similar colors 
and specific star form ation rates (SFR/M^), s uggesting similar star formation histories and 
stellar populations. In lBasu-Zych et al.l (120071 ). we found evidence for recent and contin uous 
star formation activi ty, similar to the results found for z= 3 LBGs by lErb et al.l (120061 ) and 
Shapley et al.l (120011 ). We argue that LBAs are most comparable to LBGs. 



The FUV< —18.25 sample by iHeinis et al.l (120071 ) is a less extreme sample (extending 
to fainter magnit udes) of UV-selected gal a xies compar e d to UVLGs. The other st ar- forming 
sampl es (BLGs, iMcCracken et al.l J2008h . ICoil et al.l J^OoJ, ICoil et al.l J2008h . iGiUi et al. 
(120071 )) are selected based on colo r and magnitude, a nd span a larger range of SFRs and 
stellar masses. As demonstrated in lHeinis et al.l (120071 ) for the local universe, the correlation 
strength appears to decrease with increasing UV luminosity. Yet, blue galaxies appear to 
have similar or lower yq than the GALEX sample, including, the UVLG sample. 



Kauffmann et al.l (120041 ) finds that around a critical stellar mass of ~ 10 M©, galax- 



ies separate into two classes: massive, red, with low SFRs versus less massive, blue, with 
high S FRs. The typical U VLG stellar mass is ~ W^-^Uq (LBAs have lO^" <M,/M0 < 
10^'^ '') (jHoopes et al.l 120071 ). While extremely star-forming and blue, these galaxies have 
masses placing them on the boundary between these two classes- suggesting, that they may 
be in the process of transitioning. This in-between stage is not well-understood as yet, but 
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mergers and environmental processes may drive the transition in this regime. 



The more massive, infrared counterpart to UVLGs are Luminous Infrared Galaxies 
(LIRGs) and the more extreme Ultra-Luminous Infrared Galaxies (ULIRGs) - dusty star- 
forming galaxies with significant emission cor ning from central Ikpc. T hey experience sig- 
nificant merger events in the local Universe ( Sanders fc Mirabel Il996l ). As yet, little re- 
search has been completed on s tudying the environments of local ULIRGs. According to 
Zauderer. Veilleux. fc Yed (120071 ). mostly z< 0.3 ULIRGs occupy the field environment, with 
few cases hav ing cluster densities of Abell richness class or 1. Studying LIRG pair frac- 



tions at z~ 1. [Bridge et al.l (120071 ) find that ~ 27% of the IR luminosity density arises from 
close pairs, while morphologically classified mergers contribute to ~ 34%- thereby, ~ 61% 
of the IR luminosity density can be attributed to early interactions (as evidenced by close 
pairs) and mergers (morphologically determined). While UVLGs are not dusty, it seems that 
UVLGs and LIRGs/ULIRGs share environmental conditions, possibly contributing to their 
high SFRs. 



4.2. Effects at Large scales: Importance of Halo mass 



Conroy &: Wechslerl (120081 ) investigate the relationship between star formation, stellar 
mass and halo mass. Assuming that a tight relationship between stellar mass and halo 
mass exists, based on observational "abundance matching", they theoretically derive how 
mass assembly and star formation history evolve with time. Based on averaged results 
(they caution that exceptions to their results are completely plausible), IConroy &: Wechsler 
( 120081 ) determine that the relationship between star formation rate (SFR) and halo mass is 
Gaussian-shaped, whose mean and normalization parameters evolve slightly with redshift. 
However the broadness of the SFR peak, spanning 11. 9< log(Mhaio/M0) <13.1 for z~0, 
suggests that a distinct characteristic halo mass leading to higher SFRs does not exist. 
Rather galaxy properties likely shift gradually across the halo mass scale. 



As discussed earlier, iKeres et al.l ( 120051 ) suggest two regimes of gas fueling and star for- 
mation efficiency, depending on halo mass: hot-mode accretion for massive galaxies and halos 
masses > 10^^'^ M©- which relate to lower SFRs, and cold-mode accretion in less massive 
galaxies, where gas is drawn along f ilamentary structure and a ble t o feed star formatio n more 
efficiently. Both of the results from lConroy fc Wechslerl (120081 ) and lKeres et al.l (120051 ) imply 
that the range of allowed halo masses for high SFRs could be broad- with the possibility 
that there is some mixing of both hot- and cold- accretion modes for halos with masses, 
Mhalo ~ W^MQ. 
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While the study by IConroy &: Wechslerl (120081 ) has been very thorough for z<2, they 



do aU ow that the situation might differ at higher redshifts. As illustrated by iHeinis et al. 



(120071 ). at high redshifts observations show a strong correlation between rest-frame UV- 



luminosity and correlation length, or halo mass (lAdelberger et al.ll2005l : lArnouts et al.ll2002 



Foucaud et al.ll2003l : iGiavalisco fc Dickinsonll200ll ). Yet, their results show that the nearly 
linear relation between UV- lu minosity, or SFR, and halo mass seems to level off, or possibly 
even reverse at low redshifts (IHeinis et al.l 120071 ) . 



Although we have studied a rare sample of galaxies, their close relation to well-studied 
high redshift LBGs allo ws us to make a direct comparison about the halo mass condition 
for LBG-like properties. lAdelberger et al.l (120051 ) find that LBG halo masses are ~ lO^^M©. 
The UVLGs have halo masses that are consistent with these values. While the errors are 
admittedly too large to be conclusive, our result hints that the halo masses in extreme 
star- forming galaxies, such as LBGs and UVLGs, are consistent with masses ~ lO^^M©. 



4.3. Small Scale Clustering: Influence of mergers 



At large scales, we measured halos masses of UVLGs, in an attempt to connect star 
formation properties with large scale processes; at small scales another impor tant environ 



mental quantity is measured- the likelihood for mergers. It has been shown by iBerrier et al. 



(120061 ) that halos do not necessarily merge in the same way that galaxies merge. At small 
scales, the more relevant quantity is the Halo Occupation Distribution (HOD), or the way 
in which galaxies occupy their halos. This information constrains galaxy formation models 
at the < 100 kpc scales. 

Studying their small scale clustering, we find that UVLGs and LBAs tend to have high 
pair fractions and merger rates with the ALL sample. HST images of the LBAs offer similar 
clues, displaying companions or ongoing interactio ns (such as tidal features and plumes) 
in the h i gh res olution images (jOverzier et al.ll2008l ). Consistent with the results found by 



Li et al.l (120081 ). the most extreme star-forming ga laxies are prefer entially found in close 
pairs. Studying LBGs at high redshifts (z~ 4 — 5), iLee et al.l (120061 ) also find evidence for 
small-scale clustering. They interpret that LBGs likely have fainter companions occupying 
the same halos, and find that the HOD model fits their data well. 



Using merger histories to connect z= 3 LBGs to present-day massive galaxies, IConselice 

(120061 ) argue that M^, > 10^° Mq galaxies undergo ~ 4 major mergers, with most of the 
interactions having occured before z=1.5. Furthermore, they show that z=0.5 galaxy with 
stellar mass ~ 10^° Mq likely originated as a 10^ Mq galaxy, whose stellar mass has grown 
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both by mergers and increased star formation activity. While the number of mergers plateaus 
for 10^° Mq galaxies by z~2, less massive galaxies seem to continue merging until z=0.5. 
UVLGs have stellar masses 10^°-^ — IO^^-^Mq, and LBAs have lower masses (10^-^ — 10^°-^ 
Mq). If U VLGs foll o w the same merger histories and evolution as more general galaxies 
studied by IConselicd ( l2006l ). the present UVLG and LBA masses suggest that their z= 3 
predecessors had stellar masses between 10^ — 10^ Mq and are probably still actively merging 
at z< 1. 

Could mergers be the main triggers for the increased star formation in these low dust, 
high sta r -form ing galaxies? Our results, and more generally, the conclusion reached by 
Li et al.l (120081 ) find that interactions and mergers do lead to enhanced star formation. 
As discussed in th e Intro duction, this concept is not new - having theoretical roots from 
Toomre fc Toomrd (jl972l ). mergers efficiently funnel gas into central regions of galaxies, pro- 
viding fodder fo r star forma tion. This theory is supported by several observational studies 
(see Reviews bv iKeell and lstruckl Jl999h l 



We can also reverse the query: do all mergers and interactions create UVLGs or LBGs? 
Mergers and interactions do enhance star formation, but these processes are also seen in 
LIRGs and ULIRGs - which are dustier systems, apparently different from UVLGs and 
LBAs. The overlap between these IR star-forming populations and UV-star-forming galaxies 
has not been studied in detail. We speculate that both of these systems may be related, as 
different phases of the merger process- with various unknown factors such as orientation, 
timescales, and kinematics adding complexity to the evolution and final merger state; or, 
while some fraction of mergers result in UVL Gs, some othe r fracti on result in dustier LIRGs 
or a separate class of galaxies, altogether. [Bridge et al.l (120071 ) has found that a larger 
fraction of LIRGs appear in the merger phase than in the pre-merger, close pair phase. So 
it is plausible that UVLGs relate to the earliest interaction phase, exhibiting high SFRs, 
low dust and distorted features. The highest restframe-UV surface brightnesses, observed 
in LBAs, corresponds to galaxies found to have a higher incidence of pairs. Unfortunately, 
further investigation about their environment is beyond the scope of this statistical study 
because of their rarity. Extrapolating from our data on local LBG analogs, it seems that 
close pairs are one relevant ingredient to their condition. However, it is unlikely that every 
merger will result in a low dust, high surface brightness, high SFR galaxy - other (unknown 
as yet) constraints must also be placed in order to produce this particular breed of galaxy. 
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5. Conclusions 

Using ttie COSMOS data to compare angular cross-correlation functions of UVLGs 
with three other samples: ALL, BLGs, and LRGs, we determine both large and small scale 
environments of UVLGs. In our analysis, we assume the following: 1.) The correlation 
function does not evolve with redshift. The full correlation function includes a term from 
evolution (See Eqns. [2] and [9]). We assume that e = 7 — 3, which desc ribes non-evolu t ion in 



the correlation function. This assumption is justified by the study by iMeneux et al.l (120061 ) 
- finding that ro evolves by < 20% over z=0.2— 1.2 range. 2.) The correlation function can 
be characterized b y a power law at large scales, and a separate power law at small scales (see 



Ross et al.l (120071 )). 3.) The size of the redshift bins do not affect correlation results. Bins 
that are too wide might make real-space correlations appear weaker by projecting across 
an excessively wide redshift range. However, we assume that dz=0.4 is not too wide, and 
check this assumption by comparing to more narrow dz=0.2 bins- finding that there is good 
agreement between the wide and narrow cases. 

We find that UVLGs inhabit halos with masses exceeding 1O^^M0. This result is further 
constrained by the additional observation that UVLGs and LRGs (which typically trace the 
most massive halos) do not often form pairs. Rather, UVLGs form close pairs with higher 
probability compared to the ALL sample forming pairs with themselves; and UVLGs pair 
with BLGs with the same probability as BLG-BLG pairs. Therefore we conclude that 
mergers are relevant environmental effects responsible for triggering high levels of SFR, 
observed in UVLGs and the high redshift LBGs. We are not able to definitively measure 
the halo mass of this sample given our large errors. However, we explore the possibility 
that halo mass may play some role in the UVLG or LBG condition- since both examples of 
UV-selected, high SFR galaxies have halo masses consistent with ~ IO^^Mq. 

[JGALEX (Galaxy Evolution Explorer) is a NASA Small Explorer, launched in April 
2003. We gratefully acknowledge NASA's support for construction, operation, and science 
analysis for the GALEX mission, developed in cooperation with the Centre National d'Etudes 
Spatiales of France and the Korean Ministry of Science and Technology. The HST COS- 
MOS program was supported through NASA grant HST-GO-09822. More information on the 
COSMOS survey is available at |http: / /www.astro. caltech.edu/cosmos We thank Michael 



Blanton for access to the IDL kcorrect (version 4.1.4) analysis package. This work has 
greatly benefitted from the careful comments and suggestions made by the anonymous ref- 
eree. A.R.B. gratefully recognizes Ian McGreer and Andrei Mesinger for their contributions 
to this analysis, and David Hogg for insightful discussions. 
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Table 1. Sample Numbers 



Sample 


z 




N 


All 


0.4 - 


0.8 


60872 


All 


0.8 - 


1.2 


107473 


All 


0.4 - 


1.2 


168345 


BLG 


0.4 - 


0.8 


7678 


BLG 


0.8 - 


1.2 


25100 


BLG 


0.4 - 


1.2 


32778 


LRG 


0.4 - 


0.8 


11175 


LRG 


0.8 - 


1.2 


15861 


LRG 


0.4 - 


1.2 


27036 


UVLG 


0.4 - 


0.8 


752 


UVLG 


0.8 - 


1.2 


2737 


UVLG 


0.4 - 


1.2 


3489 


LBA 


0.4 - 


0.8 


50 


LBA 


0.8 - 


1.2 


232 


LBA 


0.4 - 


1.2 


282 
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Table 2. Fits for Luminosity Function Used to Derive n(z) 



kjdiiipie 


z 


Q. 


A/T* 


/A* 


Reference 








yi\n) nidg 


j yiy} Vordl iVipC ) 






All 


0.3 


-1.3 


-21.10 


32.26 


Faber et al. 


(^007) 


All 


0.5 


-1.3 


-21.30 


33.32 


Faber et al. 


f2007l 


All 


0.7 


-1.3 


-21.62 


32.16 


Faber et al. 


Q007) 


All 


0.9 


-1.3 


-21.35 


41.26 


Faber et al. 


r2007l 


All 


1.1 


-1.3 


-21.48 


34.72 


Faber et al. 


f2007) 


BLG 


0.3 


-1.3 


-20.84 


24.26 


Faber et al. 


(^007) 


BLG 


0.5 


-1.3 


-21.20 


23.20 


Faber et al. 


f2007l 


BLG 


0.7 


-1.3 


-21.40 


27.27 


Faber et al. 


(20071 


BLG 


0.9 


-1.3 


-21.20 


37.32 


Faber et al. 


(20071 


BLG 


1.1 


-1.3 


-21.35 


29.19 


Faber et al. 


(20071 


LRG 


0.3 


-0.5 


-20.73 


21.91 


Faber et al. 


(20071 


LRG 


0.5 


-0.5 


-20.87 


19.97 


Faber et al. 


(20071 


LRG 


0.7 


-0.5 


-21.20 


17.75 


Faber et al. 


(20071 


LRG 


0.9 


-0.5 


-21.28 


11.89 


Faber et al. 


(20071 


LRG 


1.1 


-0.5 


-21.68 


5.32 


Faber et al. 


(20071 


UVLG 


0.3 


-1.5 


-19.53 


6.00 


Arnouts et al. f2005) 


UVLG 


0.6 


-1.5 


-19.63 


10.30 


Arnouts et al. f2005l 


UVLG 


1.0 


-1.5 


-19.89 


8.00 


Arnouts et al. f2005l 



'^Mg values from iFaber et al.l ( 200711 have be e n con verted from Vega mag- 
nitudes to AB magnitudes using IWillmer et al.l (j2006|). 
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Table 3. Large-scale Auto-correlation Results 



Sample 


z 




(10-=^(°)'') 


5 


ro 

(h-^Mpc) 


ALL 
ALL 


0.4 - 
0.8 - 


0.8 
1.2 


5.0 ±0.9 
22.0 +^8 


64 +°-°^ 

•J-"^ -0.04 
r, oo +0.04 
^■'->^ -0.03 


3.3 ±0.4 
3.8 ±0.4 


BLG 
BLG 


0.4 - 
0.8 - 


0.8 
1.2 


16.6 

50.0 ™ 


0.44 

n 00 +0-08 

"J-^^ -0.05 


3.9 t\l 
4.4 


LRG 
LRG 


0.4 - 
0.8 - 


0.8 
1.2 


2.7 

22.3 t^i 


0.83 
0.43 


2.9 1^:^ 
4.8 1?:^ 


UVLG 


0.4 - 


1.2 


QO Q +541.0 
_2i_i 


n 09 +0.22 
'^■^^ -0.29 


4.8 li^/ 



Table 4. Small scale UVLG Cross-correlation Results 



Samnle 


z 




(10-3(°)^) 


5 


(h-^Mpc) 


(Mpc-3) 




.'^merge 

(Gyr-^) 


UVLG X ALL 

i;VLG X ALL 


0.4 - 

0.8 - 


0.8 

L2 


n Dfi +°-^'^ 

^■^'^ -0.01 

U.UO _QQ2 


1.5 ±0.2 

-0.1 


-■-•"5 -0.5 

1 o +0.6 
-■-■■J -0.1 


0.084 ±0.006 

0.159 ±0.027 


0.47 

1.18 


+0.12 

-0.08 
+0.31 
-0.28 


0.109 t'o'oll 
n 27fi +n-079 

^' -0.064 


UVLG X BLG 

UVLG X BLG 


0.4 - 
0.8 - 


0.8 

L2 


-0.07 

n 24 +°-^^ 


1 c; +0.5 
^■^ -0.4 

1.2 ±0.2 


1 4 +3-8 
-■-■^ -1.0 

1 q +1-8 

-0.8 


0.011 ±0.001 
0.039 ±0.012 


0.07 
0.14 


+0.13 
-0.03 
+0.08 
-0.06 


U.UIO ^0.006 

n 932 +° °^^ 

U.UOZ -0.013 


UVLG X LRG 
UVLG X LRG 


0.4 - 
0.8 - 


0.8 
1.2 


0.21 

0.16 trd 


1 9 +0.8 
J--^ -0.9 
T 9 +0.3 
^■^ -0.3 


1.6 tlf 

1 S +3-0 


0.015 ±0.001 
0.023 ±0.010 


0.03 
0.11 


+0.07 
-0.01 
+0.09 
-0.06 


0.007 IHi^ 
0.025 t'of,l 


LBA X ALL 
LBA X ALL 


0.4 - 
0.8 - 


0.8 
1.2 


n 99 +10.41 
V.Z.Z. _Q 2l 

01 +°-^^ 

'J-'J-'- -0.01 


1 q +0.6 
-'-■'J -0.7 

-•-•O -0.4 


1 q +21.0 

^■^ -1.5 

1 2 +^-^ 

-•-•^ -0.9 


0.084 ±0.020 
0.159 ±0.060 


0.65 
2.61 


+1.04 
-0.,32 
+52.64 
-25.78 


n 1 p-o +0.242 
U.iOO _o.074 

0.608 ™' 



Table 5. Small scale Auto-correlation Results 



Sample 


z 




A^.e 

(10-3(°)^) 


S 


To* 

(h-^Mpc) 


n 

(Mpc-3) 


^pair 


Rmerge 

(Gyr-1) 


ALL 


0.4 - 


0.8 


0.56 t'^il 


1 OS +°0'' 

■'-•'J'J -0.05 


^••J -0.3 


0.084 ±0.015 


0.14 ±0.03 


0.033 IHS? 


ALL 


0.8 - 


L2 


0.30 ±0.05 


1 19 -O "'' 


1.9 ±0.2 


0.159 ±0.027 


0.26 ±0.05 


060 +"-"^2 

U.UUU _o.012 


BLG 


0.4 - 


0.8 


^■^^ -0.35 


1.14 ±0.24 


-■-•y -1.0 


0.011 ±0.002 


0.03 ±0.01 


007 +" ""^ 

U.UU/ _o.003 


BLG 


0.8 - 


1.2 


0.34 t'oit 


-, +0.10 
-0.11 


1 q +0.9 
-L.y -0.5 


0.039 ±0.010 


0.07 ±0.03 


Q Q17 +0.006 
-0.006 


LRG 


0.4 - 


0.8 


0.38 tri 


-1 90 +0.15 
-0.14 


2 +^-^ 


0.015 ±0.002 


0.06 ±0.02 


0.015 t^ot 


LRG 


0.8 - 


1.2 


0.18 tVo^ 


1 41 +0-08 
J-.^-L -0.09 




0.023 ±0.004 


0.24 ±0.06 


0.056 t'o'olt 
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Volume-limited Sample Selection 
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Fig. 1. — Sample selection for the three comparison samples, drawn from the I-band selected 
(I<25.5) COSMOS catalog: ALL (upper panels), BLG (middle left), and LRG (middle right) 
and the color-magnitude distributions of UVLGs (lower left) and LBAs (lower right), selected 
from the COSMOS+ GALEX (I<25.5 and NUV<23) combined catalog. The upper left panel 
illustrates the luminosity-redshift distribution for the full data, with the black box defining a 
volume-limited sample (ALL). The plot on the upper right shows how this sample of galaxies, 
ALL, occupies the color-magnitude diagram. The lower plots show color magnitude diagrams 
for the other samples. We define BLGs to have — <1.8 and Mg <-20.0 (shown on 
middle left) and LRGs to have - M, >2.0 and M, <-20.5 (middle right). The lowest left 
panel shows the UVLG distribution in purple with contour levels representing the number 
of galaxies in color-magnitide bins corresponding to (3, 6, 12, 24) galaxies; similarly, the 
lowest right panel shows the same for LBAs with contour levels corresponding to (1.5, 3, 6) 
galaxies. 
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Fig. 2. — Redshift distribution of all the separate samples - ALL (top left), BLG (top 
middle), LRG (top right), UVLG (bottom left), and LBA (bottom right). For each sample, 
we divide the full redshift range into dz=0.4 bins: 0.4—0.8, and 0.8—1.2. Black dash- 
dotted lines show predic ted number densities calculated from luminosity functions, given by 
Faber et all J2007h fand lArnouts et all J2005h . for UVLGs). 
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Fig. 3. — Spatial sky distributions-0.4<z<0.8 (left), and 0.8<z<1.2(right) for samples from 
top to bottom: ALL (top), BLG, LRG, UVLG, and LB A (bottom). The distribution is 
smoothed to a scale of 4.8'. The spatial distributions are displayed over the mask (black 
denotes regions where no reliable observations exist)- the allowed regions of the masks, 
appearing in these plots as dark violet regions underlying the colored points, are used for 
populating the random catalogs. Colors in this figure progress from deep violet to red with 
increasing density. 
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Fig. 4. — We compare the auto-correlation results for the narrow (dz=0.2) to wide (dz=0.4) 
redshift bins, by displaying the ratios of w{6), at top, and Wp(rp), at bottom. For simplicity, 
we have chosen to display only one representative redshift range (0.4—0.8) for the auto- 
correlation of the ALL (left), BLG (middle) and LRG (right) samples. The ratios of the 
0.4<z<0.6 data to the wider 0.5<z<0.9 are shown as black points, and the ratios of 0.6—0.8 
to 0.5—0.9 appear as orange points. While the narrow bin data appears a factor of 1.5—2 
higher than that of the wide bins in the w(^) plots (top), we show that converting this into 
Wp(rp) units shows little discrepancy between the separate narrow redshift ranges and the 
wider redshift- (the dashed red line marks where the narrow and wide bins give equivalent 
results. In most cases, the correlation function derived using wide z bins corresponds to the 
average correlation function of the narrower bins. Therefore, selecting a dz=0.4 redshift bin 
will not systematically affect the correlation strength. 
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Fig. 5. — Fits to the angular correlation function, for both auto-corrlelations (top) and 
cross-correlations with UVLGs (bottom) for the ALL (left), BLG (middle) and LRG (right) 
samples. The shaded regions within the dashed lines provide errors caused by cosmic vari- 
ance. The dotted vertical line delineates the physical scale of 200 kpc, the critical distance to 
separate large scale from small scale. The blue dashed lines show the large scale fit, including 
the IC correction to the data, while the red long-dashed lines mark the best fit to the small 
scales. Again, we demonstrate our fits only for one representative redshift sample: 0.8—1.2. 
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Fig. 6. — Fits to the angular correlation function, for auto-corrlelations of UVLGs. The 
shaded regions within the dashed lines provide errors caused by cosmic variance. The dotted 
vertical line delineates the physical scale of 200 kpc, the critical distance to separate large 
scale from small scale. The blue dashed line shows the large scale fit, including the IC 
correction to the data, and green solid line marks the best fit to the entire range of scales. 
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Fig. 7. — Comparing correlation functions for our samples, for the 0.8—1.2 redshift data. 
The shaded regions within the dashed lines provide errors caused by cosmic variance. The 
topmost panel compares auto-correlation of the ALL sample (black) with BLG sample (open 
navy circle) and LRG sample (open red circle). The second panel displays the ALL auto- 
correlation (black) with the 'ALL x UVLGs' cross-correlation (purple). The third panel, 
similarly, shows the BLG auto-correlation (blue) in comparison with the 'BLG x UVLGs' 
cross-correlation (purple). Finally, in the bottom panel the auto-correlation of 'LRGs' is 
shown in red, with the cross-correlation of 'LRGs' with UVLGs shown in purple. While on 
small scales, the auto-correlations of the ALL and BLG samples are indistinguishable from 
their cross correlations with the UVLGs, the LRGs appear to cluster with themselves more 
strongly than with UVLGs. 
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- UVLG correlation length, determined from the autocorrelation of 0.4<z<1.2 
We compare with correlation length data for star-forming ga laxies from other 



Fig. 8.- 
UVLGs. 

studie s. We also include high redshift results for LBGs, BX, and BM (see lAdelberger et al. 



(120051 ) for their description) with low-redshift GALEX results to compare UVLGs with 
other UV-selected samples (all shown in magenta). The dark g reen shaded region show s the 
evolution of a 10^^-^— 10^^-^ Mq halo with redshift (identical to Adelberger et al. ( 20051 )). 
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Fig. 9. — Evolution of pair fractions-Top panel: pair fraction of ALL sample with itself 
(cyan) compared to ALL with UVLGs (black). UVLGs have an increased likelihood of 
interacting with ALL galaxies than the ALL sample interacting with itself. Middle panel: 
BLGs interact with themselves (cyan) similarly to UVLGs interacting with BLGs (black). 
Finally, LRGs are marginally more likely found in pairs with other LRGs (cyan) compared 
to being found near UVLGs. 
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Fig. 10. — Comparing auto-correlation results for the three comparison samples- ALL (black 



points), BLGs (blue points), and LRGs (red points) with other published results. iLin et al. 
(I2OO4J ) (green shad ed region) s tudie s the evolution of pair fraction using galaxies with 
- 22 <Mr < -20: IBcU et all J2006h (magenta star) studies a sample with Mb < -20; 



Conselice et al.l (120031 ) (orange square) finds merger fraction, assumed to be half the pair 
fraction, for Mb < —18. These magnitude limited samples, with no color selection, are most 
comparable to our ALL sample. From Eqn. [18] we see that the pair fraction depends on 
number density, which depends strongly on limiting magnitude. In order to compare our 
pair fraction with other studies, we extrapolate the pair fraction of the ALL sample for 
Mb < —20, changing only the number density, and assuming tha t the correlation lengths is 
the same as our ALL sample (See discussion in iBell et al.l (120061 )). These points are shown 
as open circles, and appear to agree well with the other studies. 



